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Q . Abstract 

^S| ' We use recent data that show a narrow peak around y/s ~ 2.37 GeV in the pn — > dir^TT^ cross section, 

5^ . with about double strength at the peak than in the analogous pn — >■ dn^Tr^ reaction, and, assuming that it is 

C^' due to the excitation of a dibaryon resonance, we evaluate the cross section for the pn — > pmr^'!r~ reaction, 

■^^ . with the final pn unbound but with the same quantum numbers as the dcuteron. We use accurate techniques 

to determine the final state interaction in the case of the pn forming a dcuteron or a positive energy state, 

^vj , which allow us to get the pn — )> pmr^n^ cross section with pn in / = and 5 = 1, that turns out to be quite 

close or saturates the experimental pn — >■ pmr'^Tr" total cross section around ^/s ~ 2.37 GcV, depending on 

the angular momentum assumed. We then parametrize a background with different methods, and the sum of 

the resonant and background contributions is fitted to present data. The resulting cross section exceeds the 

experimental results in the region of the resonant peak, showing a problem in the dibaryon hypothesis. Yet, 

y ' in view of the dispersion of present experimental data, and the scarce information around y/s = 2.37 GeV, 

pH . a call is made for precise measurements of the pn — > pmr^ir' reaction around this energy, to further clarify 

this issue. 
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I. INTRODUCTION 

The pn — )■ dvr'^Tr'' reaction has shown an intriguing feature since the cross section exhibits a 
very narrow peak around ^/s = 2.37 GeV of about 70 MeV [1, 2]. The invariant mass distribution 
also shows a preference for the two pions having an invariant mass close to two pion masses. The 
quantum numbers of the reaction demand that the tt^tt^ state is in isospin / = which in turn also 
calls for Ltt = of the pair of pions to allow them to go together and have an invariant mass closer 
to the threshold of the two pions. The narrowness of the signal has prompted the authors of [1, 2] 
to claim that it is most probably due to the formation of a dibaryon resonance. The search for 
dibaryons has been a recurring subject, so far filled only with negative results. The confirmation 
of a dibaryon as being responsible for the reaction would undeniably be an important discovery. 

The first problem one faces in the dibaryon interpretation is that the pn — )■ dir^ir^ reaction is 
done as a fusion process. The pn scattering is actually done using a deuteron target and arguments 
are given in [1, 2] to show that the proton in the deuteron is acting as a spectator. Hence, the 
reaction proceeds via the breakup of the target deuteron followed by the recombination of the 
deuteron from the neutron of the target and the proton of the beam, according to [1, 2]. One 
might think that it is easier to make the scattering on the original deuteron that is not broken 
and remains in the final state, and actually this is an unavoidable part of the mechanisms of the 
actual reaction, although the cuts made to demand that the proton of the initial deuteron acts 
as a spectator should make that mechanism sub dominant. Yet, it is unclear what its unavoidable 
interference with the dominant mechanism could be. 

One can look for other possible reasons for the narrow peak. Assuming the fusion reaction to 
be responsible for it, if one has a long range t-channel mechanism for the production (imagine 
for instance A A production mediated by pion exchange), the fusion reaction involves the deuteron 
wave function in momentum space in a way that could magnify certain kinematics. The fact is that 
although possible explanations could be given related to the the way the deuteron wave function 
enters the fusion reaction, no theoretical study along these lines has been done. The hypothesis of 
the dibaryon resonance then stands without a contradiction so far. 

One step forward to show possible contradictions was given in Ref. [3]. In this paper the 
authors take the experimental cross section of the pn — >■ dTr^ir^ reaction, assuming it to be formed 
by a resonance in the entrance channel, and relate it to the cross section of the pn — t- pmr^Ti^ 
reaction with the positive energy pn system in the final state having the same quantum numbers 
as the deuteron, / = 0, S = 1. Using approximate techniques to take into account the final state 
interaction in the case of d or pn formation, the authors can determine the resonant cross section for 
pn — )• pmr^TT^ reaction, having a strength at the peak of about 0.4 mb.^ To this cross section, one 
should add the background contribution from many other quantum numbers in the final state and 
compare with the total experimental cross section for the pn — )■ pmr^-K^ reaction. Unfortunately 
there are no data for this reaction and the authors urged the experiment to be performed to clarify 
the issue. 

In [3] the consistency with the inelasticity in the NN cross section was also investigated and 
it was found that the inelastic NN cross section based only on the resonant mechanism exceeded 
the experimental cross section of the ^D^ partial wave obtained from the SAID analysis of [4], and 
was barely below the sum of the two partial waves ^D^ +^ G3, which are possible with J = 3^ in 
the entrance channel. In between, the normalization of the data used for the pn — >■ dn^ir'^ reaction 
has been reduced from 0.4 mb at the peak to 0.27 mb in [5] and now the estimated NN inelastic 



The authors of [3] study the case where the angular momentum between the pair of pions and the deuteron is 
L = or L — 2. Since L = 2 is the case favored in [1, 2] we refer to the numbers of [3] for L — 2. 



cross section from the pn — )■ pmr^TT^ reaction with pn having the deuteron quantum numbers alone 
would be well below the experiment.^ 

In between, some important experimental information has appeared from the measurement of 
the pn —7- d-K^TT~ cross section around the peak of the pn — )■ dir^ii^ one. Since the vr+vr" state 
has also / = component, it was certainly a puzzle that the narrow peak seen in the pn — )■ d-K^ir^ 
cross section was not seen in the pn — t- dTr^vr". Yet, the high precision data for this latter reaction 
measured in [5] also show a clear and narrow peak exactly at the same position as the pn — t- dir^iT^ 
reaction, with a strength at the peak about twice as big as for pn — )• dir^iT^ , as demanded by isospin 
symmetry. 

With this latter information it becomes most advisable to conduct the same test as in [3], 
relating the pn — )■ (ivr+vr^ and pn — )■ pnvr+Tr" cross sections, because now the predictions for the 
latter reaction can be compared with the experimental data [6-9]. On the other hand, one can 
also use a more elaborate model to account for final state interaction than used in [3], which was 
based on the use of the pn scattering amplitude in the deuteron channel as being dominated by 
the deuteron pole. In the present paper we shall study the final state interaction using techniques 
developed in the chiral unitary approach [10-14]. Using unitarity in coupled channels and the N/ D 
method it was found that it was justified to make an on shell factorization of the potential and 
the t-matrix in the Bethe Salpeter equation [12, 15] which renders this equation an algebraic one 
for each partial wave. Furthermore, we shall also use results from [16-18], where the link between 
scattering amplitudes of the chiral unitary approach and wave functions in coordinate space is 
made and a meaning is found for the couplings of a bound state or a resonance to the interacting 
particles. Using these theoretical techniques and the recent data for the pn — )• d7r~^TT~ reaction 
from [5], we shall see that the predicted cross section for the pn — >■ pmT~^7r~ reaction with pn with 
the deuteron quantum numbers, assuming it to be due to the dibaryon formation, nearly equals 
the experimental total cross section for this reaction at the energy where the pn — ^ dir^-K^ cross 
section peaks. By adding a necessary background, which is obtained by a fit to the pn — )■ pmT~^7r~ 
data adding the resonance contribution and a background parametrized in different ways, the total 
cross section obtained exceeds the experimental one around the energy of the resonant peak. While 
this poses a problem to the dibaryon hypothesis, we also show that the data are scarce around 
^/s = 2.37 GeV, and we urge that precise experiments around this region are performed to clarify 
this issue. 



II. FORMALISM 

In Refs. [1, 2] the dibaryon resonance {R) is assumed to have J = 3~^ with an orbital momentum 
L = 2 between the vr+vr" and the deuteron. In principle, L = could also be acommodated but 
analysis of angular momentum in [1, 2] favours L = 2. As in Ref. [3], we shall study both cases. 

In Fig. 1 we show the Feynman diagrams that lead to the formation of the deuteron (d) in the 
pn — 7- d'K'^TT" reaction or to the final state interaction (FSI) of the pn pair in the pn — )■ pmr~^TT~ 
reaction. In the first case (upper diagram), two nucleons are produced from a bare R — >■ pmT~^ 7r~ 
vertex, which interact among themselves to form the deuteron. In the second case (lower diagrams), 
we have the bare R — )■ pmr'^Tr" , corresponding to no interaction of the pn pair, and the further 
interaction of this pair. Note that the interaction of the pn pair is needed to form the deuteron. 

Let Vp be the bare vertex for R — >• pmT~^TT~ . We shall assume, as in [3], that this corresponds 
to a short range process, as one expects from a compact object of a resonance. Following the idea 
of [1, 2] that this object is related to a double A state, the two A's would be bound by 93 MeV, 
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FIG. 1: Fcynman diagrams leading to the formation of the deutcron in the pn — >■ dir^Tr reaction (top) and 
to scattering oi pn with positive energy (final state interaction) in the pn — >■ pmr~^T:~ reaction (bottom). 

which certainly makes the system rather compact. The transition matrix, tp, from pn to 7r'^7r~d 
is then given by:^ 



t? = VpG{Md)g , 



(1) 



where g is the coupling of the deuteron to pn {g"^ is the residue of the pn scattering matrix at the 
deuteron pole) and G is the pn propagator or loop function of the p and n propagators evaluated 
at the deuteron mass, M^. The G function is given, in the case of two nucleons, by [11]: 
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with W the CM energy of the d or the final pn system, and M the average nucleon mass. Upon 
Cauchy integration of the q^ variable, G{W) can be written as: 



r(w\= f "''^ [My 1 

^^ i (27r)3 U(g)/ W-2Eiq)+ie' 



(3) 



which is conveniently regularized with a cut off ^max by means of ^(^max" l^l)- The G function that 
appears here is the same that appears in the Bethe-Salpeter equation (BSE) (Lippmann-Schwinger 
with relativistic propagators), which in its on shell factorized form [12, 15] reads: 



t = v + vGt , t 



(-•-«)" 



(4) 



where v is now the potential and t the scattering matrix. Note that in the ordinary Lippmann- 
Schwinger equation v and t of the second term on the right hand side of the first form of Eq. (4) 



We neglect the small D-wave component of the deuteron. 



would go inside the integral of Eq. (3) with its half off shell form. The on shell factorization, where 
one neglects the left hand cut of the dispersion relation (which can be easily acommodated with a 
suitable small change in the cut off), allows to factorize v and t with their on shell value outside 
the integral, leading to Eq. (4). One can actually perform a fit in the pn data and determine v and 
the cut off providing a good fit to experiment, but we shall not need it here. 
The factorized form of the BSE is rather useful since it allows us to write: 

t = v(l + Gt) , 1 + Gt = - . (5) 

V 

Let us call t of Eq. (4) tpn^pn in what follows for clarity. If we go to the case of unbound pn 
production in Fig. 1 (lower panel), the transition i-matrix will be given by: 

i^") =Vp + VpG{W)tpn,pn = Vp (1 + G{W)tp^,pn) = Vp^^ . (6) 

Now it is useful to make use of the fact that the tpri,pn scattering matrix has a pole at the deuteron 
mass, where we have: 

tpn,vn{W) = ^^_ ^^ , W^Ma. (7) 

Now, if we look at the expression of tpn,pn in Eq. (4), we can see that if there is a pole at VF = M^ 
then: 

v~^ = G{Md) . (8) 

Furthermore, upon use of L'Hopital rule on the last form of Eq. (4), we can see that: 

1 



^^^^^ -dW 



(9) 

W=Ma 



assuming that v is energy independent, which is a quite good approximation to obtain phase shifts 
by means of Eq. (4). 

By using Eq. (8) we can now go back to Eq. (6) and write: 

4""^ = VpG{Md)tpn,pn ■ (10) 

Now, comparing Eqs. (10) and (1), we find that: 
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Hence, in order to relate the cross sections for pnir'^TT" and dn^n^ production all we need is to 
know g and tpn,pn- 

The value of g can be easily obtained from Eq. (9), and one can show that, in the limit of small 
binding, it is independent of the value of (?max [16]. In fact, using the result obtained in Ref. [16] 
adapted to our Field Theory normalization, we find: 

g' = ^^,j = V2i^=y/MB, (12) 

where n is the reduced mass of the pn system, ^ = M/2. Above, B is the binding energy of 
the deuteron. This result is actually well known since it was first discussed in Ref. [19] in the 



context of the deuteron being a simple bound state of a proton and a neutron. This result has also 
been widely used to investigate if certain states qualify as dynamically generated states from an 
interaction potential or have a very different nature [20, 21]. 

In our normalization, the tpn^pn scattering matrix is related to the one of Quantum Mechanics 
by: 
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(13) 



with k the momentum of the particles in the CM frame. In order to be as model independent as 
possible, we take tpn,pn from experiment using the pn phase shift of [4]. As we can see, the ratio 

tp /tp has dimensions MeV~ , as it should be, since in the evaluation of the cross sections 

there is an extra particle in the pn — >■ pmr~^7:~ reaction and hence an extra SpM/E{p) integration. 
The cross section for the pn — >■ (ivr+Tr" reaction is given by: 
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with s the Mandelstam variable for the initial pn system. Upon integration over the variables of 
the two pions, Eq. (14) gives: 
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where we have denoted by Mjnv the two-pion invariant mass, which is related to the deuteron 
energy by: 



M- 



s + MJ- 2^sEd . 



In Eq. (14), p is the pion momentum in the rest frame of the two pions, 

2Afinv 

being m the pion mass, and pd is the deuteron momentum in the global CM frame, 
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(18) 



In these equations, A(x, y, z) is the Kahlen or triangle function. Analogously, the cross section for 
the pn — >■ pmr~^Tr" reaction is given by: 
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which upon integration of the two pion momenta gives: 
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Here, Mjnv is the pn system invariant mass, 
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2M^ + 2EpEn - 2\pp\\pn\ cos ( 
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with pp {Ep) and p„ {En) referring the the momentum (energy) of the nucleons in the final pn 
system. Minv is the two pion system invariant mass, now given by: 
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(22) 



In order to obtain the ratio of cross sections we need 
and (13), is given by: 

2 



tp /tp , which, according to Eqs. (11) 
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with k given by: 



k 



2Minv 



(24) 



For the case L = we shall make use of Eqs. (15) and (20). In the case L = 2, however, we 
must take into account a factor q^ inside the integrals, where q is the relative momentum of the 
two-pion system with respect to the deuteron or the pn system. This amounts to put the factor 



for a^'^^ in the integrand of Eq. (15), and {\pp\ + \pn\ + 2 \pp\ \pn\ cos 6*) in the integrand of 



\Pd\ 

Eq. (20) for a^P"). 

It is interesting to compare our approach with that of Ref . [3] . The results of this latter work are 
obtained by making the approximation of the scattering length for the amplitude tpn^pn of Eq. (11), 
and demanding that it has a pole at the deuteron mass. Namely, A; cot 5 in Eq. (23) is substituted 
by -1/a with l/a = 7 of Eq. (12). 

III. RESULTS 



1.5 



0.5 



L = 

L = 2 

Dakhno et al. I — B — I 

Brunt et al. I — ■ — I 

Tsuboyama et al. I — A — I 



2.15 



2.3 2.35 

7^ (GcV) 




FIG. 2: Cross sections for the pn — > pmr^Tr reaction with pn having the deuteron quantum numbers. 
Dashed line: prediction for L = 2; solid line: prediction for L = 0. Data are taken from Refs. [6-9]. 



Once the ratio o^'P'^\^) / a'^'''\^) is calculated from Eqs. (15) and (20), u'^^") can be calculated 
from 0"^'^' i^fs)- We parameterize the experimental cross section for pn — )■ d7r"^7r~ in terms of a Breit- 
Wigner with mass Mr = 2.37 GeV and width T = 70 MeV, with a peak a{^ = Mr) = 0.44 mb.^ 
The cross sections that we obtain for pn — )■ pmT'^'K~ with the final pn having the deuteron quantum 
numbers is shown in Fig. 2 for L = and L = 2, together with the experimental data for the total 
cross section of pn — )■ pnir^Tr^ [6, 7]. As we can see, for the case L = the predicted cross section 
is barely below the experimental point at ^/s ~ 2.37 GeV. For L = 2 the cross section is about 
25% below the datum. 

It is interesting to compare these results with those found in Ref. [3] for the pn — )■ pmr^n^ 
reaction. The results are qualitatively similar, the cross section for L = being larger than that 
for L = 2. The ratio of cross sections a^P"-) /a^'^' is 0.9 in Ref. [3] for L = 2, while in our case it is 
1.2 for the charged pions reactions, and it would be the same for the neutral pion reactions. The 
qualitative agreement is quite good in view of the simplifications done in Ref. [3], as quoted above. 
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FIG. 3: Cross sections for pn — >■ pmr^iT^ when a background is added to the resonant cross section. The 
solid line includes a background determined through phase space, while the dashed line includes a polynomial 
background. Data are taken from Refs. [6-9]. The left panel shows the whole range 2.15 — 2.55 GeV, whereas 
the right panel contains the energy range yG = Mr ± F. 

The results of Fig. 2 indicate that if we add some background for the pn — )■ pmT~^TT~ reaction 
from all other contributions where the pn does not have the deuteron quantum numbers, or from the 
standard non-resonant mechanisms of two-pion production present in theoretical models [22, 23],^ 
one would get sizeable cross sections that might exceed the experimental cross sections at the 
peak of the resonant contribution. In order to show this, we have made some estimation of the 
background in two ways. In one case, the background is taken from phase space. This is easily 
obtained, up to a constant factor, from Eq. (20) setting tp ftp to unity in the integrand. In the 
second case, we parameterize the background by a polynomial in the variable P = \/s — 2M — 2m, 
that is: 



<7BG 



CiP + C2P2 + CsP^ 



(25) 



where q are arbitrary constants. We perform a best fit to the data suming the background and 
the resonant cross sections. The results so obtained are shown in Fig. 3, only for the case L = 2. 



* The value at the peak in Ref. [5] is 0.5 mb, but to this one must subtract the background from 1 = 1 making also 

small corrections for different phase space for charged and neutral pions. 
^ The model of Ref. [23] does not sum amplitudes for the different diagrams but cross sections. Still, it provides fair 

integrated cross sections in the different channels. 



We can see that in both cases (phase space and polynomial backgrounds) the cross section around 
the resonant peak exceeds the experimental cross section. Yet, we also observe that there is some 
dispersion of the data and the information around the relevant region of y/s = 2.37 GeV is scarce. 
In view of this, we can only encourage measurements of the pn — t- pmr~^Tr~ cross section with good 
resolution in the range y^ = 2.30-2.45 GeV in order to clarify the situation and show wether there 
is or not a strong resonant peak. 
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FIG. 4: Predicted cross sections for the "pn — >■ [dlpn)'K^'K^ reactions. The blue solid line with circles 
represents a Breit-Wigner parameterization of the experimental cross section for pn — >■ di^Pi^P of Ref. [5], 
with a = 0.27 mb at the peak. Our prediction for the pn -^ pn-n^-n:^ reaction (having the final pn the 
dcutcron quantum numbers) is shown with the red dashed line. The black solid line results from the sum of 
both contributions. The upper dot-dashed lines represent the SAID analysis [4] for the NN inelastic cross 
section for the ^D-^ and the ^D-^ +'^ Gj, partial waves, with an isospin factor 1/6. 



We also perform the same test as done in Ref. [3] to compare the inelastic cross section for 
NN scattering with / = with the cross section obtained for the pn — )■ pmr~^7r~ reaction. From 
isospin considerations the / = 0, J^ = 3"^ contribution from A^A^ — )■ NNtttt would be three times 
the cross section of pn — t- pmr'^Tr" or six times that of pn — t- pnir^iT^, to which we should add the 
contribution from pn — t- cItt'^tt' or pn — )■ dn^ir^, respectively. In view of this, in Fig. 4, we plot the 
cross section for pn — t- dTr^n^ and pn — )■ pmr^-K^ and their sum as compared with one sixth of the 
inelastic cross sections for A^A^ in / = 0, J = 3~^. Two possible initial partial waves can contribute 
in this case, the ^D^ and the ^G^ waves. We show the inelastic cross sections corresponding to 
the SAID analysis [4]. We see that the cross section evaluated from the two pion production is 
somewhat below the contribution from the ^D^ and the sum of the ^D^ and ^G^ waves. It should 
be noted that, in spite of getting a larger ratio for a^P"-' /a^ ' here than in [3] (1.2 versus 0.9) we 
still get a smaller A^A^ — )■ NNtttt cross section. This is because the 0.4 mb cross section used in [3] 
at the peak was changed to 0.27 mb in [5]. Thus, the argument used in [3] to put constraints on 
the resonance hypothesis is weakened when the new data are used. 



IV. CONCLUSIONS 



We have carried out an exercise, taking into account the final state interaction of a pn pair with 
the deuteron quantum numbers to produce a deuteron or a pn pair with positive energy, in order to 
relate the cross sections of the pn — >• d'7r~^7r~ and pn — >■ pmr^7r~ reactions with the final pn pair with 
the deuteron quantum numbers. The test was done assuming that the clear narrow peak seen in the 
pn —7- dir'^TT" reaction was due to the formation of a dibaryon resonance in the pn entrance channel. 



as suggested in the experimental papers of [1, 2, 5]. We found that the cross section for this latter 
reaction was about 1.2 times bigger than for the pn — >■ d-Tr+vr" reaction, assuming that one has L = 2 
angular momentum between the pn and tttt pairs, as claimed in the experiments. We compared 
this cross section with the total experimental one for the pn — )■ pn-Tr+vr" reaction and found that it 
was quite close to the experiment at the peak of the resonance. We then added a background from 
unavoidable standard nonresonant channels, which summed to the resonant contribution was fitted 
to present pn — )■ pnvr'^vr" data, and found that the total cross section exceeds the experimental 
one around ^/s = 2.37 GeV, where one has the peak of the pn — )■ dvr+vr" cross section. We could 
admit uncertainties of about 20% in the size of the pn — )■ pmr~^TT~ cross section (as also admitted 
in Ref. [3]), but the fit with the resonant plus background contributions leads to results basically 
the same as in Fig. 3. Hence, based on this scarce experimental information, this excess puts 
a problem to the hypothesis of the dibaryon resonance as being responsible for the narrow peak 
observed in the pn — >■ dTT^Tr^ and pii — )■ dTT~^n~ reactions. In view of this, we can recommend two 
lines of research to further clarify this problem. The first one is to work on theoretical models for 
A^A^ — >■ NNirn, and concretely for the NN — )■ dnir, reactions. For some mechanisms which involve 
a light particle (a pion) exchange in the t-channel, the consideration of the deuteron wave function 
is bound to constraint some distributions in the phase space that can repercute also in the cross 
section as a function of the energy. Some steps in this direction were done in [24] in the study of 
the pn — )• dn^TT^ reaction using the model of [22], but at much lower energies than those discussed 
here. The second line of suggested research is the measurement of the pn -^ pmr~^7r~ cross section 
with good resolution around the region of ^/s = 2.37 GeV, where the peak of the cross sections 
in the deuteron production processes is seen. We have clearly shown that a large resonant cross 
section should develop as a consequence of the assumption of a resonant formation in the entrance 
channel. So far, the scarce data around this energy do not let us see the shape as a function of 
energy. A combination of efforts in those or other directions to clarify why such a narrow peak 
appears in the deuteron fusion reactions are much needed to see if the resonance hypothesis stands. 
We should note that many of the data points are very old, and present experimental facilities at 
COSY, HADES and other can provide measurements with far better accuracy. In view of this, we 
urge that experimental efforts are devoted to this task which, combined with the findings of the 
present paper, should help clarify this interesting problem. 
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